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We consider the behavior of an Ising ferromagnet obeying the 
Glauber dynamics under the influence of a fast switching, random ex- 
ternal field. After introducing a general formalism for describing such 
systems, we consider here the mean-field theory. A novel type of first 
order phase transition related to spontaneous symmetry breaking and 
dynamic freezing is found. The non-equilibrium stationary state has a 
complex structure, which changes as a function of parameters from a 
singular-continuous distribution with Euclidean or fractal support to 
an absolutely continuous one. These transitions are reflected in both 
finite size effects and sample-to-sample fluctuations. 
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I. INTRODUCTION 

The last decade has seen many advancements in the theory of dynamic systems, 
both from a mathematically rigorous, as from a physically oriented heuristic point 
of view. Most of these results have been obtained for systems with few degrees 
of freedom, attempts at handling nonequilibrium stationary states of sjfstems with 
macroscopically many degrees of freedom have been made only recentlytTH. In this 
paper we propose a general theoretical framework for strongly, randomly driven 
statistical physical systems and apply it to the Ising model in a random, dichotomic 
driving external field. 

The resulting dynamics has many qualitative similarities with earlier work on one 
dimensional Ising chains in a binary random fieldP. However, while the 'dynamics' 
defined ind is a one dimensional map generated by the iteration of non-commuting 
2x2 transfer matrices, in the present case the map 'lives' in 2 N dimensions. Nev- 
ertheless, the main mechanism leading to chaotic behavior and strange attractors 
is in both cases related to the competition between two (or more) fixed points (or 
limit cycles). In this respect, the randomly driven Ising model (RDIM) has many 
intriguing aspects which - somewhat unexpectedly - can be handled both analyt- 
ically fftd, numerically with methods developed earlier for the random field Ising 
modelErM 

The Ising ferromagnet in a time-dependent sinusoidally oscillating field has re- 
ceived recently a lot of attention, from both a theoretical and experimental point 
of view. On the theoretical side, Rao, Krishnamurthy, and PanditOJ have presented 
a large TV-expansion of the cubic 0(N) model in three dimensions and calculated 
the critical exponents related to the area of the hysteresis loop. The underlying dy- 
namic phase transition has been then studied within both mean-fieldEJ and Monte 
Carlo simulationsoO. The theory presented in this paper is a generalization of 
these ideas for the case when the external field is subject to a chaotic dynamics 
and/or is a random variable. 
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Besides the theoretical interest in describing such systems, we believe that many 
of our predictions can be tested with recently developed experimental techniques. 
Dynamic magnetization measurements have been recently performed io. ultrathin 
Au(lll)/Cu(0001)/Au(lll) sandwiches or epitaxial Co/Au(lll) glmsBil Sim- 
ilarly, hysteresis measurements on the ultrathin film Co/Au(001)E2l indicate that 
below T c these systems undergo a dynamic phase transition belonging to the Ising- 
universality class. More relevant to our theory, the time evolution of magnetization 
clusters can be optically recorded. The typical relaxation times range from minutes 
to a few seconds with increasing field amplitudes^. This relatively slow relaxation 
rate allows for a simple experimental realization of the randomly driven external 
field. 

Ultrathin films are potential candidates for magneto-optical storage devices and 
our approach might be relevant especially in this respect. At well chosen control 
parameters the stationary magnetization distribution of the RDIM displays several 
well separated peaks. Thus, when driven appropriately, such materials can store 
locally more than the two values typical for an equilibrium ferromagnetic system. 

This first paper is organized as following: the basic assumptions and the gen- 
eral theoretical formalism are introduced in Section 2. The mean-field theory is 
presented in full detail in Section 3. In this approximation the paramagnetic- 
ferromagnetic stationary phase transition becomes first order. The phase boundary 
is obtained analytically. The average magnetization is non-analytic (jumps) at small 
driving field values and therefore the usual mean-field approach fails. Nevertheless, 
the phase transition is related to a spontaneous symmetry breaking and a pitchfork 
bifurcation of the magnetization distribution. On the other hand, the analytical 
nature of the stationary magnetization distribution changes from singular contin- 
uous with euclidean or fractal support to absolutely continuous along analytically 
computed boundaries. Such changes are directly connected to finite size effects of 
the free energy and the multifractal spectrum of the magnetization measure. Close 
to the para-ferro phase boundary the small magnetization region close to m ~ be- 
comes a repellor, inducing ultracritical slowing down related to type-I intermittency. 
A short summary of the main mean-field features is presented in Section 4. 

A subsequent paper considers the randomly driven Ising model in one and two 
dimensions. Although the one-dimensional case cannot be fully solved, many inter- 
esting exact results can be derived. We find a line of second order phase transitions 
at T — between a disordered (driven paramagnetic) and an ordered (ferromag- 
netic) phase. Along this line the stationary critical exponents change continuously 
as a function of the driving field strength. However, the dynamic critical expo- 
nent remains unchanged, z = 2. As a function of temperature and of the driving 
field strength, the nonequilibrium stationary state might display multifractal or 
'fat' multifractal character. Hence, the generalized free energy is characterized by 
anomalous fluctuations related to the existence of a multifractal spectrum. We 
also performed Monte Carlo simulations on a square lattice. Many features of the 
mean-field dynamics are shown to survive the strong fluctuations characteristic to 
two dimensional systems. However, in contrast to the mean-field approach, the two- 
dimensional model displays also an interesting spatial structure related to droplet 
dynamics. Comparisons between mean- field theory and two-dimensional results are 
systematically presented, including some preliminary results for hysteresis. 



II. GENERAL FORMALISM 
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A. The Master Equation 



Consider a spin system fx — (s%, S2, ■ ■ ■ , sat) in contact with a thermal bath 
(phonons). Let us denote by r sv i n -fu v the characteristic time of the spin-phonon 
interactions, and by T sys the slowest relaxation mode of the spin system. For 
Tspin-fUp •C T sys the system is in local thermal equilibrium. 

Our basic assumption is that the time evolution of the spin system can be de- 
scribed by a time-dependent joint probability distribution, P({si};t). From a dy- 
namic point of view, one can consider this distribution as an average over all mi- 
crostate initial conditions satisfying given macroscopic constraints. The rigorous 
definition of this basic assumption- is related to the existence of Markov-partitions 
in (hyperbolic) dynamic systemsEm and goes beyond the scope of this article. 

The distribution \P(t)) obeys the Master Equation: 

\P(t)) = -£ m \P(t)), (1) 

where the operator C-B(t) describes the outflow (inflow) probability from (into) state 
{s^ into (from) other states, as prescribed by the dynamic rules. The dependence 



on the field is made explicit because in general 



^ for B 7^ B' . Assuming 



a constant field B = Bq and ordering the eigenvalues of Cb as Aq = < | A x | < 
••• < |A 2 n_ 1 |, one has T sys = |Ai| _1 . The ket vector \P(i)) can be expressed in the 
spin-configuration basis as P({si};t). Alternately, it can be parameterized as 

1 N 1 ( 2N \ 

i—l i^j \ a — 1 jGct J 

where the 7r a (£) are the average values of all possible 2 N products of spins 

* a (t) = (l[s j ) = Y / P({si};t)l[s j (3) 

j£a s 

Hence, \P(t)) can be expressed as a 2 Ar -dimensional normalized vector in the 
space of all possible spin configurations, or in the space of spin-products as 
7?(t) = ((si), . . . , (sjv), (S1S2), . . . )• Note than wheiu. expressing the kinetic Ising 
model Liouville operator in terms of Pauli matricesEHI this orthogonal transforma- 
tion corresponds to the exchange o~ z <-> o~ x . 

In general, the Liouville operator C is not symmetric but can be expanded in a 
biorthogonal basis formed by its right, \r n ), and left, (l n \, eigenvectors ((l n \r, n ) ~ 

^m,n) ■ 

2 N -1 

t= Yu i r «) A «^i ( 4 ) 

n=0 

It is worth noting that since e~ c is a stochastic operator, we have (lo\C = 0(Zo|, 
where (Zo| = (1)1,..., 1) in the spin-configuration basis. The scalar product 
(^o| r o) = Z delivers the equilibrium (stationary) canonical partition function. 

The time-dependent field B(t) is usually a deterministic one-dimensional map, 
e.g. a harmonically oscillating field. However, if the deterministic map is chaotic, 
the field becomes a random variable. In what follows we assume that the exter- 
nal field B is a random variable sampled identically and independently from the 
symmetric distribution p(B) = p(—B). 



3 



Let tb be the average sampling time of the field distribution. As long as tb 3> r sys 
the spin system has enough time to relax to global thermal equilibrium. This is 
the case of equilibrium statistical mechanics and normal fluctuations. The situation 
can be very different if the field can switch abruptly {r~^ itch ~ /iflB is large) and 
i~sy S 3> tb- The system does not have enough time to relax to equilibrium and 
the stationary state is determined by the external random field. As discussed later, 
such a situation is experimentally realizable. 

Now let us assume that the field is sampled from p(B) at time intervals of length 

TB, 

OO 

B(t) = p B Bp{B) 9(* ~ nr B )Q ((n + l)r B - t) (5) 

n=0 

We can integrate Eq. (|]) exactly from t n -\ = (n — 1)tb to t for t n -i < t < t n : 

\P({ Si };t)) = e-^-^-^lPdsi};^)) (6) 

where B(t n -i) is the field instance sampled at time t n -\. For t = lim e ^o(t n — e) 
one obtains 

\P({s t }; t n )) = e-^n-O-B | F({Si}; tn _ i)} (?) 

Since in our case the low lying eigenvalues A Q of the Liouville operator C satisfy 
AqTb f, by expanding the exponential in first order one obtains the discrete 
'coarse grained' Master Equation 

\P({ Sl };t n ))-\P({ Sl };t n ^)) _ t \ P(U \.t V (ffi 

- - L B(t n - 1 )\P{{Si},t n -i)) (a) 

tb 

which describes correctly only the long-time behavior of Eq. ([I]). Short term effects 
due to the larger eigenvalues of C have already relaxed at the time scale rg. This 
approximated form of the Master Equation is used in all further developments. 

One can regard Equation (0) as defining the discrete dynamics governing the 
probability distribution \P({si\;t n )) . 



B. The driving field distribution 



As already mentioned, the external field might be distributed according to the 
invariant measure of some chaotic one-dimensional deterministic map. In other 
applications, the field can be Poisson- or Gauss-distributed. In what follows we 
restrict ourselves to the binary distribution 

p(B) = ±5(B-B ) + ±6(B + B ) (9) 

Many of our results can easily be generalized to arbitrary continuous distributions. 
Other results, in particular those concerning the stationary state phase diagram and 
the critical behavior, depend strongly on the discrete character of the choice (^). 
Eqs. (||) and (^) map our problem into an iterated function system (IFS)Eix As 
long as we deal with a finite system of spins, the mathematical results (including the 
collage theorem) developed by Demko and Barnsley and subsequent work on IFS 
apply to randomly driven spin models as well. However, from a statistical physical 
point of view, the interesting things happen after taking the thermodynamic limit. 
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C. The invariant measure 



As usual in the theory of dynamic systems, one caB-.ask what is the invariant 
measure induced by the dynamics (||), see, for example,Eil. Let V s (jf) denote the in- 
variant density related to the dynamics Eq. ||. It satisfies the Chapman-Kolmogorov 
equation: 

V s (tt) = J dn' P s (tt') J dB p(B) 5(n - C - &BTB vr') = JCV s (n) (10) 

where K. denotes the Frobenius-Perron (FP) operator. Physically, V s (tv) describes 
the nonequilibrium stationary state induced by the Master Equation dynamics. 
Note that we used above the spin-product basis. An orthonormal transformation of 
the basis will lead to a different but equivalent FP-operator, the transformation's 
Jacobian is unity. Again, this is true only for finite systems. 

If one is interested in the stationary expectation value of some spin observable 
A({si}) one must perform two averages, the 'thermal' (...) and the 'dynamic' average 



A({ Si }) = [(A)] (11) 

where the 'dynamic' average [...] is taken over V s and p{B). 

The 'thermodynamics' of such driven systems can be computed from the gen- 
eralized free energy. This is related to the largest Lyapunov exponent A of the 
dynamics as — f3J- = A. Consider a long dynamic trajectory consisting of Ttb 
sampling points. The Lyapunov exponent is defined as 

A = lim — !— l n T I ! e -^B( T )TB e -£. B (,T-i)TB ... e -c B (i)T B \ ( 12 ) 
t^oo tbT I 1 

where B{n) = B(t n ) is distributed according to Eq. ([)]). The same result can be 
obtained by iterating some general^ initial unity vector \po) as 
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|p 2 ) - — c- Cb ^ Tb \ Pi ), a 2 = VHpaj 

For large n the vectors \p n ) will be distributed according to V s and up to 0(-^) 
corrections the Lyapunov exponent can be expressed as 

A = J dV s (n) J dB p{B)± In || C -^7r|| (13) 

For a constant field p(B) — 8{B — B ) and the stationary distribution is V s {tt) — 
(5(7? — TT eq ), where n eq are the Boltzmann-distribution averaged spin-products. There- 
fore, ||e~ £B o rB 7? e9 || = ((Zo|r )) 2 = Z 2 , where we have used that n eq is the right 
eigenvector of Cb„ with eigenvalue 0. We recover the usual definition of free energy 
by multiplying A with -fc^T. 

This generalized free energy might display anomalous fluctuations related to the 
multifractal spectrum of the stationary distribution V Sl as will be shown later for 
the mean field theory. 



1 This vector should not fall into any invariant subspace of both C±b operators 
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D. Dynamical properties 



In order to consider the dynamical properties of randomly driven systems one 
has to solve - in full analogy to the theory of one-dimensional maps - the right 
eigenvalue problem of the Frobenius-Perron operator: 

The largest magnitude eigenvalue is one, sq = 1. The right eigenvector TZq is 
nodeless and real: it corresponds to the stationary state, TZq = V s . For m > the 
eigenvectors R m (x) satisfy 

J dn R m (7r) = (15) 

In addition to the right eigenvectors TZ m , the operator K, might have also a set of 
null functions of different orders. They correspond to zero eigenvalues of /C 9 , where 
q is an integer (the order). These eigenfunctions, however, do not contribute to the 
relaxation of the initial probability distribution towards V s . 

In analogy to the usual transfer matrix theory, the relaxation of the probability 
distribution and of the time dependent correlation functions are determined for 
asymptotically long times by the second largest eigenvalue s\ and its eigenfunction 
Tlx. 



Although at this stage the formalism looks rather involved, it is a straightforward 
extension of the methods developed for low-dimensional dynamic systems. We con- 
sider next the simplest possible example, an Ising model in a random binary external 
field, Eq. ([)]) . In this case, many interesting stationary and dynamic properties can 
be obtained analytically, or with numerical methods not more complex than those 
used for one dimensional maps. 

III. MEAN FIELD APPROXIMATION 

A. The Mean-Field Map 

Consider an Ising model defined on an A^-dimensional simplex, such that all spins 
are nearest neighbors: 

E = -j^^SiSj - n B B^Si (16) 

i=jLj i 

J is normalized so that the energy is additive and [ib is the Bohr magneton. 

Let jli := (si, . . . , —Si, . . . , sjv)- We may describe the Liouville operator C with 
the transition rate w{fli\(l) in the Glauber formed 

- - Sitanh(- s i + H )\ ( 17 ) 

where (3 = 1/kgT, K = f3J, H = P^bB and a sets the time constant. Applying 
Eq. (||) one obtains after performing the thermodynamic limit N — > oo: 

m(t + 1) = tanh(Km(t) + H(t)), (18) 
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Time is measured in units of tb- The field distribution Eq. (|9|) leads to the one- 
dimensional map 

!taxth(Km(t) + Hq) with probability \ 
(19) 
tanh(ifm(t) — Hq) with probability i 

Note that in the thermodynamic limit the moments of the magnetization do 
not couple with higher order correlation functions and the methods worked out 
previously for the one-dimensional random-field Ising chain can thus be applied 
directly. 

Since in the stationary state, Eq. (|To|), [m k (t + 1)] = [m k (t)}, using Eq. (18) and 



simple algebraic manipulations we obtain that the fc-th moment of the stationary 
magnetization is given by 



\m k ] 



v 



1 + vh 



fc = l,2,... (20) 



where v = tanh(Km) and h = tanh(TJ). 

At high temperature the system is in the disordered, paramagnetic phase, in 
which case all odd moments of the magnetization vanish. Assuming that the free 
energy is analytic in [m], the critical temperature is obtained by expanding [m] in 
first order in 0([ft. 2 ]): 

[m]~K{l-hl)[m] (21) 

and neglecting [m q ], q = 3 and higher odd moments (which should scale as 0([h 2q ])). 
In the usual mean-field scenario [m] = in the paramagnetic phase and the coeffi- 
cient vanishes at the transition point to the ferromagnetic phase: 

rr(II) _ 1 l T , 1+mt (00\ 

Hc 2 ln T^t ' (22) 



where — ±y K K 1 for K > 1. Using a simple geometric argument we will show 
below that this analyticity assumption fails and the phase transition is actually first 
order. For the second moment one obtains 

1™ 2 1 = i-^(i-lg + 3 ^) < 23 » 

where we have omitted 0([m 4 ]) and higher even moments. The pole of this expres- 
sion is also related to the phase transition, which is discussed below. Third and 
fourth order expansions of [m] and [m 2 ] read 



and 

\m 2 ] ~ [ //, 



h\K^{2 - llhl + 30feg - 15/tg) 

~ 3 - K 4 {3 - QQhl + 212/ig - imh% + 105/ig) 

i _ v*h - Ah 2 _i_ , K 6 (2 - \7hl + 30ht - 15fog)(6fo 2 - 16feg + lOfrjj) \ " 
1 4Al o + ^oJ+ 3 _ ^4(3 _ 60^ + 212ftg-260fcg + 105fc§) / 



(25) 



respectively. A high order expansion of the moments along these lines can be easily 
obtained using algebraic manipulations programs but will not be presented here. 
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B. The stationary phase diagram 



In principle, there are at least two different mechanisms for a phase transitions in 
the stationary state described by V s (to) . The first one corresponds to spontaneous 
symmetry breaking leading to a continuous phase transition. In this scenario the 
stationary distribution, which at high temperature is a function of the even magneti- 
zation moments only, V s (m) = V s (—m), becomes degenerate at certain parameter 
values {K,Hq} and the odd subspace, P (m) = —P (—m), contributes as well. 
More precisely, V s (m) = TZq is always a nodeless even function of to. The field 
symmetry is spontaneously broken when the largest odd-subspace eigenvalue of the 
Frobenius- Perron operator Eq. (jlj), sj. — * 1. Therefore, the largest eigenvalue is 
degenerate and the corresponding eigenvector is an arbitrary linear combination of 
V s (m) and lZi(m) = —TZi(—m), leading to a non-vanishing order parameter. Close 
to but above the transition point the relaxation time of the stationary distribution 
diverges as r _1 ~ 1 — si. We find no evidence for such a mechanism, at least not in 
mean-field approximation. Instead, the phase transition is related to a bifurcation 
of the stationary magnetization distribution. 

Consider the map Eq. (19) at high temperature, a situation shown in Fig. [l|. The 



arrows indicate the direction of the flow. The competition between the two stable 
fixed points leads to chaotic behavior and the displayed stationary distribution. 
To approximate the distribution, we tracked the evolution of 1000 (random) initial 
values of to subject to the map for 1000 iterations. At low temperatures (K > 1) and 
large fields one has the situation depicted in Fig. 0. Note the possible intermittent 
behavior close to m ~ 0. If we decrease Hq the map can 'pinch' tangentially the 
m(t+l) = m(t) diagonal, creating thus one new unstable fixed point. This situation 
is shown in Fig. |^. Decreasing the field even further, we have the map of Fig. [|, 
where V s has bifurcated into two stable and one unstable disjoint distributions. 
For further use let us denote by toi , to 2 , and 7713 the possible fixed points of the 
equation to = taiih(Km + Hq) in descending order. The line of the critical field 
H c can be calculated from the condition that at the new fixed point the map is 
tangential ('critical map', seec3) and leads after elementary calculations to 

1 1 - TO 1 * 

H c = -In— — + Km) (26) 

2 1 + TO T 

where to 1 ~ to 2 — TO3 . 

From Figs. |^-[| it is evident that |to| > to 1 and except for Hq = the magne- 
tization jumps at the phase transition. We believe that this feature is due to the 
discrete character of the binary p(B) distribution. Thus, the RDIM provides an 
example of a spontaneous symmetry breaking leading to a first order phase transi- 
tion. The mechanism behind this first order transition is very different from that 
of equilibrium systems and is related to a tangential bifurcation of the stationary 
distribution. The corresponding phase diagram is shown in the upper part of Fig. 

I 



C. The multifractal regime 



Some highly unusual properties of the RDIM are related to the multifractal spec- 
trum of the stationary state. Following the notation introduced ir£3, one can iden- 
tify a singular-continuous density with fractal support (SC-F) in both the param- 
agnetic and the ferromagnetic phase. When a gap opens between the upper and the 
lower branch of the map the invariant distribution has a fractal support with the ca- 
pacity dimension do < 1. The border of the (SC-F) region is given by Km\ = Hq in 
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the para- and K(m\ +VTI3) = 2Hq in the ferromagnetic phase. In the region between 
do = 1 .and doo = 1_ the distribution is .singular-continuous with Euclidean support 
(SC-E)El. Using the ideas developed mtJ, we obtain d 00 = lii K(l - m{) = |. The 
density distribution is absolutely continuous (AC) if all generalized dimensions^ 
equal one, d q = 1, (a = 0, . . . , 00). These results are graphically summarized in the 
lower part of Fig. H. 

In order to compute the generalized free energy, Eq. ([l3|), one can use that 
(lo\ r n) = Sq.u J2{ Si ] Peq({si})- Note that the left eigenvector (Z | of the Cb opera- 
tor corresponds to a sum over all spin configurations and is therefore independent of 
B. When inserting in the product within the trace of Eq. (|l^) the spectral decom- 
positions of different non-commuting operators Cb, the ground state contributions 
decouple from the higher level contributions. Hence, in mean field approximation 
the free energy is given as expected by 

-PF = N J dmP s (m)~ In 2[cosh(2A"m) + cosh(2# )] (27) 

This integral can be approximated above H c by expanding the integrand in even 
moments of the magnetization (see Eq. ( p3| ) and (25)). The fluctuations of the free 
energy depend on the (multifractal) structure of the stationary V s (to) distribution. 

Strictly speaking, Eq. ( |27j ) is the average free energy. When considering a fi- 
nite system or a long but finite dynamic trajectory, the free energy is normally 
distributed. As shown irE for the one dimensional random field Ising model, in the 
SC-F region the multifractal spectrum can be directly related ta.the second cumu- 
lant of the free energy distribution. The arguments presented irJ3 apply also to our 
case, a broad multifractal distribution leads to large free energy fluctuations. 

In the SC-F regime one can obtain additional information about finite-size free- 
energy fluctuations from the generalized dimensions (Lege^dje^tcansform of the 
multifractal spectrum) . We applied the methods developed irfj&EJo and computed 
numerically the multifractal spectrum of the stationary distribution. 

Another interesting observation is that these isolines cannot directly cross into 
the ferromagnetic region: close but above the phase transition there is no positive 
gap (see Figs. H~H). Nevertheless, in the ferromagnetic phase the magnetization 
distribution itself can be multifractal. This is shown in Fig [l| The inset shows 
the enlarged part of the map leading to a multifractal distribution for positive 
magnetization (a symmetric counterpart exists for negative magnetization). 



D. Dynamical properties 



The stationary phase transition at H c , Eq. (26), is from a physical point of view 
a dynamic freezing transition characterized by an extremely slow dynamics. As 
shown below, the relaxation of the map - and hence of all time-dependent correlation 
functions - diverges exponentially fast close to the critical field H c . Consider first 
the mean-field map close but above the critical field, as illustrated in Fig. ^. 

The iteration along the upper branch alone corresponds to type.. I .iatermittency 
and has been discussed previously in the theory of chaotic mapscJu'Eil. As usual, 
the function m' = t&nh(Km + Hq) will be approximated up to quadratic order close 

to the point (m^,_ff c ) where rrv = ±y Z is the point where the upper branch 



touches tangentially the m! = m line. Introducing the new variable x = ™ '"!■>, , 

J iv(l— m T ) 7 

one obtains 

1 9 -^0 — He , , 

x n+ i = x n + m'x„ H — (28) 

A 
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Requiring that x " + g Xn , x 2 , and H ° K Hc have the same order of magnitude implies 
that Sn and hence n must scale as 



H — H c 
K 



(29) 



which is the standard result for one-dimensional mapscj. 

However, the probability to stay on the upper branch of the map for n consecutive 
steps is exponentially small. Assume that at time t = one injects No points at 
the m = — 1 location. In order to move upwards, the points can use only the upper 
branch and must pass through the 'intermittent tunnel'. Once a trajectory flips 
to the lower branch, it is set back to the entrance of the tunnel. If a point has 
passed through the tunnel, it might eventually return to the lower part but has a 
similar chance of being trapped on the symmetric upper part. This dynamics can 
be modeled by the Markov process shown in Fig. [?]. 

By iterating the corresponding stochastic matrix (or by full induction) it is easy to 
see that the stationary probability of being at site n is given by p n = ^ . Therefore, 
assuming quasi equilibrium, the escape rate is estimated as 

N(t) = ~N(t) (30) 

where a is a constant of order 0(1) related to the probability of return after escape. 
The relaxation time r corresponding to Eq. ( |30|) diverges as 

t = — — ~ 2 a ^- H ^~ h (31) 
a 

where we have used Eq. (p9|), a is a constant. Hence, the relaxation time diverges 
exponentially fast close to the phase transition. Below iJ c , the slow dynamics is 
due to the average escape time (fractal dimension) from the central repellor. 

Another interesting dynamic phenomenon is the hysteresis of the RDIM. Here 
one adds a harmonic part to the external driving field: 

H(t) i-f H(t) + A cos Qt (32) 

The resulting hysteresis distribution is shown in Fig. ^. The evolution of 500 initial 
values was tracked for 2000 iterations. We close here the discussion of the mean-field 
(or infinite dimensional) RDIM. We expect many of the features discussed here to 
be valid in three dimensional systems and to a lesser extent in two dimensions. 



IV. SUMMARY AND DISCUSSION 



In this paper we have discussed the behavior of a spin system with short range 
interactions in a random external field coupled to the order parameter. If the distri- 
bution of the external field is discrete, the resulting dynamics is chaotic due to the 
competition between different equilibrium states of the system. We proposed a gen- 
eral formalism for calculating the stationary and dynamical properties of randomly 
driven systems and applied it to the Ising model. In the mean-field approximation 
the stationary distribution of the magnetization displays a spontaneous symmetry 
breaking phase at low fields and temperatures. The transition between the disor- 
dered and the ferromagnetic phase is first order and corresponds to a tangential 
bifurcation of the underlying map. Close to the phase transition the characteristic 
relaxation time diverges exponentially - leading to dynamic freezing. Depending on 
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the control parameters, the stationary magnetization distribution can be a normal, 
multifractal or fat-fractal distribution in both the disordered and ordered phases. 

Our interest in this problem arises mainly in connection to understanding the 
nature of open systems with many degrees of freedom. Information processing 
systems, natural or artificial, have a macroscopic number of connected elements 
subject to external stimuli changing faster than the characteristic thermal relaxation 
time. As illustrated by the simple example presented in this paper, such systems 
might develop stationary states far from equilibrium which might be many times 
more effective in dynamically storing information than simple thermal equilibrium 
states. In this respect it would be also of interest to consider other choices for the 
driving field distribution. Continuous distributions, for instance, might lead to very 
different stationary phase transitions than the one discussed here. 

In "The randomly driven Ising ferromagnet" , Part II, we discuss the RDIM in 
one and two dimensions. 
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FIG. 2. Mean field map and the stationary distribution for K = 2.0 and H /K — 0.5. 
The map shows that for strong driving fields the system remains in the paramagnetic phase 
even below the equilibrium critical temperature. 
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FIG. 3. Same as in Fig. 1 but close to the critical field value (K — 2.0 and 
Ho/K = 0.266). Two disjoint distributions are created around the stable fixed points, 
a repellor in the middle. 
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FIG. 4. Mean field map and the stationary distribution in the ferromagnetic phase 
(K = 2.0 and H /K = O.l). 
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FIG. 5. Mean field phase diagram. The upper part (Ho > 0) shows the border between the 
para- and ferromagnetic phase. In the lower part (Ho < 0) the regions denoted by SC-F, 
SC-E correspond to a singular- continuous invariant density with fractal and Euclidean 
support, respectively, while in the AC-region the density is absolutely-continuous. Note 
that the diagram is actually symmetric in Ho ■ 
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FIG. 6. Lower part of the mean field map for H > H c . Iteration along the upper branch 
is type I intermittent while the lower branch brings the iteration to the starting point in 
one step. 
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FIG. 7. The Markov process representing the passage through an intermittent tunnel. 
Each arrow indicates a transition probability of \. 
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